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Appendix A. Pyrometer Measurements with a Hole Pre-Cut in the Booster 

Cover Simulator 

To more accurately assess the response of the high-speed video and the infrared pyrometer to test 
conditions, especially when the test produces a hole through the booster cover simulator (often 
referred to as the diaphragm), a special test was perfonned. 

In this test, a hole with a diameter of 0.1 15-inch was cut in the center of the booster cover 
simulator. The booster cover simulator is a 304 SS disk 0.003 inches thick. The pre-cut hole 
diameter was chosen as a reasonable approximation of some of the holes produced in single 
NASA Standard Initiators (NSI) tests with the stainless steel (SS) V-Primer Chamber Assembly 
(PCA). Note that area visible through the sapphire window is 0.250 inches in diameter. Thus, 
the pre-cut represented about 46 percent of the visible diameter and about 2 1 percent of the total 
visible area. 

The sapphire window and perforated diaphragm were assembled in a SS V-PCA. The V-PCA 
has flow paths with the nominal 0.060-inch diameter. Assembly was performed according to the 
standard procedure used for all the tests. A single NSI was fired in the SS V-PCA to see if the 
pyrometer and high-speed video would give reasonable results. 

Figure A-l shows a graph of the temperature of the underside of the booster cover simulator as 
seen by the infrared pyrometer. Note that the first visible video indication occurs at 123 ps 
(microseconds) after the first application of current to the NSI. Other data indicated that the NSI 
end closure opened at 1 12 ps, which would make the first picture 1 1 ps after the end closure 
opened. The camera takes 20,000 frames per second, so frames are 50 ps apart. Also note that 
the upper detection limit for the pyrometer is 2,000 °C or 3,632 °F. 

The sequence of pictures shown in Figures A-2 to A-6 show light and then hot gases arriving 
first at the booster interface followed by what appear to be burning particles of zirconium 
potassium perchlorate (ZPP). The white circle in Figure A-2 is one of the 0.060-inch-diameter 
flow paths from one of the NSI chambers in the top of the V-PCA. Part of a circular arc can be 
seen along the left-hand side of the picture. This is part of the edge of the pre-cut hole in the 
booster cover simulator. 

SS melts at about 2,550 °F. The metal changes color from a darker red-brown to yellow or white 
as it approaches that temperature. The last picture is after all the ZPP has burned and the only 
light is coming from the cooling, white-hot SS. 

In the picture taken at 973 ps, the larger, irregular circle with some areas still glowing is the 
melted edge of the pre-cut hole in the diaphragm (see Figure A-6). To the left and down, a 
portion of a still larger circle can be seen. That is the bottom or outer edge of the sapphire 
window and the “booster” cavity. The smaller circle with the glowing edge is the flow passage. 
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Figure A-7 shows a graph of the NSI current and cavity pressures versus time for SS V-PCA 
with 0.1 15-inch-diameter pre-cut hole in the diaphragm. 
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Figure A-l. Booster Interface Temperature versus Time, SS V-PCA with 0.1 1 5 -Inch-Diameter 

Pre-Cut Hole in Diaphragm 


NESC Request No.: 09-00596 






© 

NASA Engineering and Safety Center 
Technical Assessment Report 

Document #: 

NESC-RP- 

09-00596 

Version: 

1.0 

Title: 

Pyrovalve Booster Interface Temperature Measurement 

Page #: 

5 of 92 




Figure A-2. At 123 /.is, One of the 0.060 NSI Flame Channels can be seen Illuminated. The 
Pyrometer is just Starting to Indicate Temperature with a Reading of 684 °F 



Figure A-3. At 1 73 /is, the Pre-Cut Hole in the Booster Cover Simulator Can Be Seen 
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Figure A-4. At 223 / is 



Figure A-5. At 373 /is 
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Figure A-6. At 973 fis, the Edge of the Now Partially Melted Pre-Cut Hole can be seen AND the Still 
Glowing Edge of One of the Flame Channels. A Larger Arc at the lower left-hand corner is Part of the 
Hole through the Retainer Nut that Holds the Sapphire Window Assembly onto the V-PCA. 
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Time Referenced To Trigger Signal (ps) 


Figure A-7. NSI Current and NSI Cavity Pressures Versus Time, SS V-PCA with 0.115-Inch- 

Diameter Pre-Cut Hole in Diaphragm 
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Appendix B. Pyrometer Noise 

During testing, a question arose as to whether interference from the firing circuitry or some other 
source might be affecting the booster interface temperature readings from the infrared pyrometer. 

Therefore, a special test was performed. In this test, two NSIs were fired with a 250-us skew. 
These NSIs were from a different lot than those used for the actual temperature measurement 
tests. This was because it was desired to perform all the temperature tests from one lot of NSIs 
and the team wanted to conserve that particular lot. For this checkout test, whichever lot was 
used would not make any difference. 

The main difference in this test was that the pyrometer lens was completely blocked. The 
pyrometer would not be able to sense any energy from the booster interface at all. If any 
disturbance was seen in the pyrometer readings, it would have to be from some other effect. 

As seen in Figures B-l and B-2, the amperages and pressures produced during the test were 
normal. However, there was no effect on the pyrometer readings. 


Therefore, it was concluded that the pyrometer readings seen during the actual test series were 
unaffected by other test conditions. 



Figure B-l. Pyrometer Readings Hold Steady at Minimum when Two NSIs Were Fired with a 250-f.is 

Skew and when the Pyrometer Lens is Blocked 
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Time Referenced To Trigger Signal (ms) 

Figure B-2. Pyrometer Readings Hold Steady at Minimum when Two NSIs were Fired with a 250-jus 

Skew and when the Pyrometer Lens is Blocked 
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Appendix C. Pressure Transducer Drop Test 

During Run 6 of Phase IIB, one of the pressure transducers failed to register any pressure 
readings at all. Troubleshooting revealed this to be a failure of the pressure transducer itself. 

As a result of this failure, a new pre-test checkout procedure for the pressure transducers was 
initiated. 

Each pressure transducer was installed in an oil-filled reservoir. A tube or pipe was installed in 
the top of the reservoir such that a weight could be positioned at a certain height and held in 
place with a mechanical pin. When the pin was removed, the weight would fall down inside the 
tube and make contact with a rod that would then push into the oil reservoir. The pressure 
transducer would then sense this pressure shock. 

Figures C-l and C-2 show typical pressure traces that are obtained from the drop test. Data were 
collected once every microsecond. The data are not suitable for calibrating the pressure 
transducers. However, the test was valuable to ensure that the pressure transducers work and the 
cable connections are good. 



Figure C-l. Typical Drop Test Results Showing Pressure Transducer Output in Volts versus Test Time 
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Figure C-2 shows the same data except that the 100 data point moving average is plotted. 
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Figure C-2. Typical Drop Test Results Showing Pressure Transducer Output in Volts Versus Test 

Time. 1 00 Data Point Moving Average. 


Figures C-3 through C-5 shows the drop tester. 
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Figure C-3. The Assembled Drop Tester. A Pressure Transducer is Installed into the Oil Reservoir at 
the Bottom. The Drop Weight is Alongside. Holes Are Drilled in the Guide Tube at Various Heights 
so that the Drop Weight can be Positioned to get the Desired Effect. 
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Figure C-4. The Drop Tester with the Guide Tube Removed to Show the Rod that Transmits the 

Impact of the Weight to the Oil Reservoir 
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Figure C-5. In this Photograph of the Oil Reservoir, the Impact Rod has been Removed to Show the 
Installed 0-Ring. This Photograph also provides a Closer View of the Pressure Transducer 

Installation. 
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Appendix D. Assessment of Area versus Temperature Indication 

During Phase IIB of the test program, a question arose over whether the larger flow passage 
diameters might be biasing the pyrometer readings by increasing the heated area within the 
pyrometer field of view (FOV). 

In other words, consider the following situation where the green circle represents the pyrometer 
FOV and the red circle is the heated area of the diaphragm. In the illustration below, the green 
circle representing the FOV is the same size. However, the red circle on the right is twice the 
diameter of the one of the left to represent the case where the flow passage being tested had 
4 times the nominal cross-sectional area. Could this “FOV effect” give a higher temperature 
reading even if the temperature is the same? 



The manufacturer’s literature gave FOV information at various distances from the target and the 
pyrometer was located a specific distance from the sapphire window of the test article to give the 
desired FOV; but, to experimentally verify that a special test was performed. 

A target heat source was positioned behind where the underside of the booster cover simulator in 
the test article would normally be. The temperature of the target was unimportant so long as it 
remained constant at a value above the minimum detectable by the pyrometer. A variable 
aperture was located in the exact plane where the underside of the test article booster cover 
simulator would be. The diameter of the aperture was incrementally reduced until the 
temperature indicated by the pyrometer started to drop off. 

Figure D-l shows the results of the four test runs that were made and how the determination of 
the FOV was detennined. 
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Aperture Area (sq. inches) 


Figure D-l. Pyrometer Signal versus Aperture Area for a Constant Temperature Target 

The diameter of the variable aperture was measured using calipers. The measured diameter was 
used to calculate the area of the aperture, which was then plotted against the pyrometer signal in 
millivolts. The objective was to look for a sharp change in slope of the curve to indicate at what 
point the aperture was starting to restrict the pyrometer’s FOV. 

The data were divided into two groups near where the change in slope occurred. A least-squares 
regression analysis was performed on both data groups to determine the best lines to represent 
the data. The intersection of those two lines gave the aperture area where the attenuation began 
to occur. From that area, the aperture diameter was calculated to be 0.22 inches. This then was 
the experimentally determined FOV. 

The measured diameter of the underside of the booster cover simulator that is visible through the 
sapphire window is 0.25 inches. The NSI flame channel diameters ranged from 0.060 to 
0.120 inches. Therefore, the pyrometer FOV was set appropriately to give the best possible 
indication of the actual temperature on the underside of the booster cover simulator. Another 
point worth noting is that the FOV was always the same from one test to another so that the 
measured temperatures would not be biased in any case. 
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Appendix E. PCA Thermal Analysis 

Heat transfer in the PCA may result from a variety of mechanisms described in detail in 
Reference 1 . In this analysis, three heat transfer mechanisms to the booster cap surface were 
studied, individually: 

1 . Convective heat transfer from hot combustion gases; 

2. Heat transfer via deposition and subsequent solidification of molten zirconia (Zr02) spray 
resulting from combustion; and 

3. Deposition and subsequent combustion of ZPP. 

Heat transfer in the PCA is complex and the analyses required simplifying assumptions to allow 
analytical study. These will be discussed further for each particular assumption. 

All thermal network models developed for this analysis were created using Thermal Desktop R 
and solved using the Systems Improved Numerical Differencing Analyzer/Fluid Integrator 
(SINDA/FLUINT B '). Both products were developed by Cullimore and Ring Technologies. 

It is recognized that there are performance differences between the SS and A1 PCA units. No 
attempt was made in this analysis to explain the differences. Rather, it is seen that both PCA 
designs require heating of the booster cap to actuate the device. The team’s focus is on 
mechanisms that affect heating to the booster cap. Compression heating was not investigated in 
detail for this study and may have a significant contribution to the overall heat transfer to the 
booster cap. Such an analysis would require a flow analysis that was beyond the scope of this 
work. 

E.1.1 Lumped-Mass Assessment of Convective Heat Transfer to the Booster Cap 

A simplified closed-fonn thermal assessment of the SS booster cap response was perfonned 
considering only convective heat transfer to the cap and the effects of its thermal mass. The 
thermal analysis network schematic is depicted in Figure E.l-1. 

T g (Gas) 
hA 

mC p (Booster Cap) 



Figure E.l-1. Simplified Lumped-Mass Network Representation 
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For this analysis, the booster cap was assumed to be a single lumped-mass heated only by 
convection. With this assumption, the overall heat transfer is given by: 

mC p y t = hA(T g -T(t )) EQ. 1 

where m is the booster cap mass, C p is the booster cap specific heat, h is the convective heat 
transfer coefficient, A is the booster cap surface area exposed to the hot gas, T (t) is the booster 
cap temperature at time, t, and T g is the gas temperature. In this analysis, the specific heat, 
convective heat transfer coefficient and gas temperature are assumed to be constant. The 
constant gas temperature assumption results in a conservative analysis since the hot gas in the 
real configuration does not have infinite capacitance and would cool rapidly as heat is transferred 
to the PC A wall and booster cap. 


9 = 


Tg-TQ : ) 

T g- T o 


where T 0 is the initial booster cap temperature. Then: 


where r is the time constant given by: 


dt t 


mC v 

T — 

hA 


EQ. 2 


EQ. 3 


EQ. 4 


Furthermore, if the top of the booster cap is represented as a disk of radius, r, and thickness, d, 
the expression for the mass becomes: 


m — pV — pnr 2 d 


EQ. 5 


where the surface area exposed to the hot combustion gases is given by: 

A — nr 2 

The expression for the time constant, r, simplifies to: 

mCp pnr 2 dC p pdCp 

hA hnr 2 h 


Rearranging EQ. 3 yields: 


dO 

e 


-dt 


T 


EQ. 6 


EQ. 7 


EQ. 8 
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and integration leads to: 


In 0 = - -t + C EQ. 9 

T 

where C is a constant of integration. But at t — 0,T — T 0 and8— 1. This implies C — 0. 
EQ. 9 becomes: 


In (Si™ 

V Tg-T 0 


) = 


th 


pdC v 


EQ. 10 


Simplifying and solving for h: 


h = 



EQ. 11 


EQ. 1 1 allows exploration of the convective heat transfer coefficient required to produce 
temperatures seen on the booster cap. By specifying a constant gas temperature, T g , and a time, 
t, at which a booster cap temperature T (t) is observed, the required h value may be determined. 
By comparing the magnitude of h required to produce the observed temperature rise with 
published literature, an assessment of the contribution of convective heating to the overall 
booster cap temperature rise may be performed. 

The h value required to produce a booster cap temperature of 866 K (1,100 °F) at a time of 
500 ps for various constant gas temperatures is presented in Figure E. 1-2. 
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Heat Transfer Coefficient, h, for Booster Cap Temperature to Reach 
866 K (1 1 00 F) in 500 ps as a Function of (Constant) Gas Temperature 
Gas Temperature Assumed Constant 



Gas Temperature (K) 


Figure E.l-2. Convective Heat Transfer Coefficient and Gas Temperature Parametric Analysis Results 


Figure E.l-2 shows the required heat transfer coefficient is far beyond that found in the 
published literature [ref. 2] for lower gas temperatures (< 2,500 K), which suggests heat transfer 
coefficients on the order of: 

1,000 W/m 2 K to 6,000 W/m 2 K (176 Btu/hr ft 2 °F to 1,056 Btu/hr ft 2 °F) 

However, for gas temperatures above 2,500 K, the range of required heat transfer coefficients is 
within the range specified in Reference 2. 


NESC Request No.: 09-00596 


© 

NASA Engineering and Safety Center 
Technical Assessment Report 

Document #: 

NESC-RP- 

09-00596 

Version: 

1.0 

Title: 

Pyrovalve Booster Interface Temperature Measurement 

Page #: 

22 of 92 


This intuitively large range required convective heat transfer coefficients in concert with the 
assumption of constant gas temperature for this simplified analysis suggests that the observed 
booster cap transient temperature response may not be explained by convective heat transfer 
alone. Therefore, additional heat transfer mechanisms are investigated. 

E.1.2 Zirconia (Zr(> 2 ) Deposition 

It was concluded in the previous section that high convective heat transfer coefficients are 
required to account for the rapid temperature rise observed in the booster cap. Another possible 
energy transfer mechanism is the exchange of heat through the deposition and solidification of 
molten ZrCF onto the PCA internal surfaces including the booster cap. 

The NSI contains 1 14 mg (2.51 x 1 0‘ 4 lbm) of ZPP. An estimated 71 percent [ref. 1, Section 
7.4. 1 .2.2.2] of this material after combustion exists in the form of ZrCF and is expected to 
condense on surfaces cooler than 3,000 K (4,940 °F). The effect of ZrCF deposition and 
solidification and the associated heat transfer was simulated in the thermal model by assuming 
the following: 

a. 71 percent of 1 14 mg (80.94 mg) (1.78 x 10‘ 4 lbm) of molten ZrCF is deposited uniformly 
onto the PCA internal walls and the top of the booster cap; 

b. The ZrCF deposition is present at the initiation time of the thermal model (t = 0 ps) and is 
in perfect contact with the A1 walls and SS booster cap top; 

c. The assumed initial temperature for the PCA is 293 K (68 °F) and 2950.3 K (485 1 °F) for 
the ZrCF which is just above the Zr02 melting temperature of 2,950 K (4,850 °F) [ref. 3]; 

d. The heat of fusion for ZrC >2 is 707 J/g (304 Btu/lbm) [ref. 4]; 

e. The internal exposed surface area of the A PCA was calculated to be 8.4 cm (1.302 in ) 
based on the computer-aided design (CAD) representation derived from Reference 5. 
Subsequently, a 0.762-cm (0.3-inch) diameter booster cap was added to close out the 
volume; and 

f. The assumed density of the ZrCF was 5.68 g/cm 3 (0.205 lbm/in 3 ) [ref. 3]. 

Assumptions (e) and (f) may be used to determine a ZrCF deposition thickness, assuming 
uniform coating of all interior surfaces, of 0.017 mm (6.657 x 10‘ 4 inches). 

A one-dimensional thermal network model was developed to study the effect of ZrCF spray onto 
the booster cap and is depicted in Figure E.l-3. 
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Figure E.l-3. Zirconia Spray Thermal Network Model 
The predicted temperature response is depicted in Figure E. 1-4. 
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Predicted BoosterCap Response to Zirconia Spray at Zirconia Solidification Temperature 
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Figure E.l-4. Zirconia Spray Thermal Analysis Results 


— Zr02 

— Zr02/Booster Cap Interface 

— BoosterCap In-Depth 1 
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— BoosterCap In-Depth 4 

— BoosterCap In-Depth 5 

— BoosterCap In-Depth 6 
BoosterCap In-Depth 7 (Bottom) 



The effect of zirconia solidification is seen in Figure E.l-4. Latent heat stored in the melted 
zirconia conducts into the booster cap causing temperatures to rise. When the zirconia has 
solidified, conduction to the booster cap continues and a temperature drop is observed as the 
flow of sensible heat. While a zirconia deposit of uniform thickness can raise the temperature of 
the bottom of the booster cap to temperatures greater than that required for booster propellant 
ignition, the timeframe required for this is well beyond what has been observed in PCA firing 
tests. Additionally, with the assumption of constant thickness deposition, there is insufficient 
energy stored in the deposited zirconia to raise the temperature of the booster cap above its melt 
temperature. Melting of a portion of the booster cap has been seen in many of the PCA firing 
tests. 

A sensitivity analysis was performed by assuming double zirconia deposition thickness. The 
same thermal model was used except that the mass of the zirconia was doubled and the 
conductance from the double-thickness zirconia node was halved to account for the increased 
thickness; no additional nodal subdivision was added. The results of the analysis are presented 
in Figure E.l-5. 
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Predicted Booster Cap Response to Zirconia Spray at Zirconia Solidification Temperature 
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Figure E.l-5. Zirconia Spray Thermal Analysis Results Assuming Double Zirconia Thickness 
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As seen in Figure E.l-5, the additional zirconia mass now contains sufficient energy to raise the 
booster cap nodes to the melt temperature. Melting progresses through the thickness until a 
temperature rise above the assumed melt temperature of 1,672 K (2,550 °F) is seen indicating 
that the melt-through was complete. While this analysis results in melting and temperatures 
sufficient to ignite the booster propellant, neither occurs within the timeframe observed during 
testing. 

Finally, it should be noted that only through-the-thickness heat transfer was modeled and 
conduction losses to the PCA structure were not modeled. In actual operation, these losses 
would reduce the energy available to heat the booster cap. 

From this simplified analysis, it can be concluded that deposition of molten zirconia onto the 
booster cap may play a significant role in raising the temperature in addition to convective heat 
transfer. However, the timeframe during which heat transfer to the booster cap occurs suggests 
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that this is not the primary mode of heat transfer to the booster cap. The sensitivity analysis 
suggests that thicker deposits of zirconia increase the propensity of melting the booster cap and 
accelerating heat transfer through the cap. 

E.1.3 ZPP Combustion in Contact with the Booster Cap 

It is estimated that up to 20 percent (22.8 mg) of the ZPP exits the NSI unburned [ref. 1]. High- 
speed photography obtained during NSI testing shows that a portion of the ZPP emerges from 
the NSI unburned as is seen in Figure E.l-6 [ref. 6]. 



Figure E.l-6. Photographic Evidence of Late ZPP Combustion 


The heat transfer associated with transport of ZPP to the top of the booster cap was investigated. 
ZPP has a heat of combustion ( H zpp ) between 5,609 and 6,070 J/g and a density (p zpp ) of 
3.23 g/cm J (90 lbm/in 3 ) [ref. 1, Section 7.4. 1.2. 2.1]. In the scenario examined, unburned ZPP is 
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expelled from the NSI and some u nkn own fraction is transported to the interior-facing surface of 
the booster cap. The purpose of this analysis was to determine whether or not this ZPP, burning 
while in contact with the booster cap, could produce sufficient heat transfer to account for the 
booster cap temperature rise, response and observed melting. Both simplified hand calculations 
and a transient thermal analysis were used to investigate this possibility. 

E. 1.3.1 Lumped-Mass Model 

In this simplified analysis, heat is distributed over the entire booster cap area exposed to the hot 
combustion gases and a transient analysis is performed to determine how much ZPP is required 
to melt the booster cap surface. The thermal model network is depicted in Figure E. 1-7. 


ZPP Heat 

O rnC p (Booster Cap) 
A Radiaton Boundary 


Kev 

A Boundary Node 
O Diffusion Node 


Figure E.l-7. Simplified ZPP Heating Model Network with Back-Side Radiation 


A parametric analysis was perfonned assuming deposition of various amounts of the initial total 
unbumed ZPP onto the booster cap surface. ZPP burning was assumed to occur over a period of 
500 ps. The results of the parametric study are presented in Table E. 1-1 . 
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Table E.l-1. Results of Parametric Analysis Investigating the Quantity of ZPP Required to Melt 

the Entire Booster Cap Surface 


Percentage of 
Unburned ZPP 

(%) 

Percentage 
of Total ZPP 

(%) 

Peak Booster 
Cap Temperature 
(F) 

Melt Status 

1.0 

0.2 

244 

None 

10.0 

2.0 

1830 

None 

11.0 

2.2 

2006 

None 

12.0 

2.4 

2182 

None 

13.0 

2.6 

2354 

None 

14.0 

2.8 

2358 

None 

15.0 

3.0 

2550 

Partial 

17.0 

3.4 

2550 

Partial 

19.0 

3.8 

2550 

Partial 

20.0 

4.0 

2621 

Complete 


With this model, a conservative upper limit for the amount of ZPP required to melt the booster 
cap was established. Note that the results are expressed as a function of percentage of unburned 
ZPP that is deposited onto the booster cap surface. If 20 percent of the ZPP from the NSI is 
unbumed [ref. 1], then 20 percent of that quantity is required to fully melt the booster cap surface 
under the assumptions in this analysis — or only 4 percent of the initial NSI ZPP load 
(0.20 x 0.20 = 0.04). This corresponds to a ZPP mass of approximately 4.6 mg (l.OxlO' 5 lbm). 

However, photographic evidence of booster caps indicates melting is a more local phenomenon 
as seen in Figure E.l-8 [ref. 7]. 
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Figure E.l-8. Phase IIB — Post-Test Images of the Booster Cover Simulators, Sapphire Windows, 
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While the results presented in Table E. 1-1 are suggestive of an upper bound to the required 
quantity of ZPP, an assessment of the temperature of local melting of metal in direct contact with 
the ZPP globule is desired. 

E.l.3.2 Local Heating from a ZPP Globule — Transient Predictions 

A refined calculation may be used to account for the physical extent of a ZPP globule by 
assuming the molten globule is a hemisphere as it sits upon the booster cap (Figure E. 1-9). 



In this instance, the mass of a ZPP globule is given by: 

m z pp — Pzpp^zpp ~ 2 Pzpp EQ. 12 

and the radius of the hemispherical ZPP globule is, then: 


r zpp = 


3 / 1.5 m 


zpp 


pn 


EQ. 13 


The resulting globule radius is used in a simplified transient thermal analysis of the booster cap 
to determine if local melting will occur, assuming that all heat generated during ZPP combustion 
is transferred to the local booster cap surface. This analysis represents an improvement in 
fidelity over the simplified lumped-mass model presented in the previous section, in that heating 
is applied only to the region in contact with the ZPP and heating of the remaining booster cap 
surface is through conduction. The added benefit of the time dependency allows comparison of 
the transient temperature profile observed during PCA testing. 

The predicted transient response to burning of a small hemispherical globule of ZPP 
(m = 6.5x1 0" 4 mg) for selected booster cap nodes is shown in Figure E.l-10. In this simulation, 
heat generated as a result of the combustion of a small amount of ZPP conducts into a 
representation of the booster cap. Heat transfer is assumed to be in the radial direction 
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(i.e., only a single node through-the-thickness). At this level of fidelity, it is clear that the 
heating is a localized effect and the temperature rises as a function of radial distance diminishes 
rapidly. This is suggestive of localized melting, an effect observed in booster cap testing 
(Figure E.l-8). Heat transfer through-the-thickness of the booster cap is investigated in the 
following section. 


Predicted BoosterCap Response to Burning ZPP, Radial Nodalization Only 
Globule Radius = 4.6x1 O' 2 mm (0.0018 in), Globule Mass = 6.5x1 O' 4 mg (1.4x1 O' 9 Ibm) 



Figure E.l-1 0. Predicted Temperature Response of Booster Cap with Localized Globule of ZPP 

(Radial Nodalization Only) 

E.l.3.3 Local Melting through the Booster Cap Thickness 

In a previous section, an estimate of the maximum amount of ZPP required to melt the entire 
booster cap was given to be approximately 4.6 mg (1.0 x 10' 5 lbm). An estimate of a lower 
bound on the amount of ZPP required to initiate local melt -through is sought. In this analytical 
development, heating resulting from ZPP burning is assumed to be transferred to the booster cap 
metal directly beneath it (i.e., a cylindrical core through-the-thickness); no radial heat transfer or 
back-side radiation is assumed. 

The mass of a hemispherical ZPP globule was showing in Eq. 12 to be: 

^zpp — ~^Pzpp (3 ™zpp ) — ^Pzpp^t^zpp EQ. 14 
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The heat liberated by burning this quantity of ZPP is given by: 

hQzpp irizppHzpp g HzpppzppTCVzpp EQ. 15 

For simplification purposes, it is assumed that the hemispherical globule of ZPP heats only the 
cylindrical core of SS directly beneath it. Then, the heat required to raise the temperature of the 
cylindrical core to the melt temperature plus the heat required to melt the material is given by: 


Expanding yields: 


A QTotal ~b A QMelt 


AQrotai ~ nr2 dPss[Cp(.T Me it T 0 ) + Hf\ 


EQ. 16 


EQ. 17 


Where C p is the SS specific heat, T Melt is the melting temperature, T 0 is the initial temperature, 
and Hf is the SS heat of fusion. 

For this simplified case, the ZPP globule will have sufficient energy content to melt the metal 
beneath it when: 


A Qzpp — AQrotai EQ. 18 

Substituting in expressions from EQs. 15 and 17 and rearranging gives the desired result: 

^ 3dpss[ c p(TMelt- T o)+ H f] in 

r > — — EQ. 19 

2 H Z ppp Z pp 

2 

When the appropriate parameters are substituted into EQ. 19, the resulting radius is 4.6x10' mm 
(1.8x10° inches). This corresponds to a ZPP mass of 6.5xl0‘ 4 mg (1.4x10° lbm). 

This result suggests that even the smallest of ZPP globules has sufficient energy to cause local 
melting and sets a lower bound to the ZPP required to melt through the booster cap. However, 
the local heating and melting assumption may not be realistic given that heat applied at one 
location on the booster cap can conduct both radially and through-the-thickness of the booster 
cap. A higher fidelity model is required. 
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E.l.3.4 Two-Dimensional Axisymmetric Model 

A two-dimensional axisymmetric thermal model was developed to investigate the effect of 
burning ZPP on a lumped-mass representation of the booster cap including localized melting 
effects with model fidelity in both the radial and through-the-thickness directions. Phase change 
for the SS was incorporated into the model with the following assumptions: 


a. Transient analysis was performed assuming transport of the unbumed ZPP to the booster 
cap surface in the form of a hemispherical globule; 

b. ZPP combustion with constant heat generation over a period of 500 ps; 

c. Heat transfer to the booster cap was assumed to be perfect; all heat generated was used to 
raise the booster cap temperature and melt the material, if applicable; 

d. No transfer from the booster cap to the PCA structure was assumed; and 

e. Radiation off the back-side is to a 293 K (68 °F) environment. 

A schematic representation of a booster cap finite difference thermal model is shown in 
Figure E. 1-1 1. 
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Figure E.l-11. Burning ZPP Thermal Network Model 
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2 3 

Heating was calculated and applied for the case where r = 4.6x10' mm (1.8x10' inches), a value 
corresponding to the radius determined by the simplified analysis in the previous section, dubbed 
Case 1 . The results of the transient analysis are presented in Figure E. 1-12. 


BoosterCap Temperatureas a Function of Depth Versus Time 
Case 1-- Burning ZPP Model, 6.5xl0‘ 4 mgZPP 



Figure E.l-12. Two-Dimensional Axisymmetric Thermal Model Results for a 4.6 xlO' 2 mm (1.8x10 
3 inches) Radius Hemispherical ZPP Globule Showing Incomplete Melt-Through 


As shown in Figure E.l-12, there is sufficient energy to raise the temperature and melt nodes on 
the upper booster cap surface. However, due to radial conduction effects plus differentiation of 
nodal temperatures through-the-thickness, there is insufficient energy to melt completely through 
the cap. Predicted temperatures on the bottom of the cap reach a maximum of approximately 
755 K (900 °F), far short of the assumed melt temperature of 1,672 K (2,550 °F). 

To determine the amount of ZPP required to melt through the entire booster cap depth, the radius 
was raised to twice the previous value r = 0.091 nun (3.6x10' inches), dubbed Case 2. Since the 
volume of the associated hemisphere of ZPP changes with r 3 , the ZPP mass for the larger radius 
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globule is 8 times that in the previous case; hence, 8 times as much energy is liberated by this 
larger globule. However, the surface area over which the heat is applied increases as r 2 . 

The math model was rerun producing unrealistically high temperatures on the booster cap 
surface — temperatures far in excess of the melting point and well above the expected ZPP 
temperature. This is due to the chosen modeling technique wherein ZPP burning is considered to 
be an applied heating rate rather than a combustion temperature with a specified conductance to 
the booster cap. However, from a heat transfer through- the -thickness perspective, the increased 
heating provides sufficient energy to melt completely through the booster cap. Temperatures 
along the bottom of the booster cap surface as a function of time and radial distance are shown in 
Figure E. 1-13. 

While not a full fidelity analysis, the result suggests that the smallest amounts of ZPP burning in 
contact with the booster cap contain sufficient energy to produce the observed melting of the SS 
cover. It is also observed in Figure E. 1-13 that melt-through is a highly localized effect. Given 
the mitigated temperature response as radial distance increases, the assumption of negligible heat 
transfer from the booster cap to the PCA structure is validated. 


Burning ZPP Temperature as a Function of BoosterCap Radial Distance vs. Time 
Case 2 -- Burning ZPP Model, 5.2x10 3 mg ZPP 
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Figure E.l-13. Two-Dimensional Axisymmetric Complete Melt-Through Thermal Model Results for 

the Booster Cap Bottom 
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The results of this refined analysis suggest that complete melt-through of the booster cap can be 

3 8 

accomplished with as little as 5.2x10' mg (1.14x10' lbm) of ZPP forming a hemispherical 
globule with r = 0.091 mm (3.6x10" inches). 

E.1.4 Conclusions 

Three heat transfer mechanisms and their effect on booster cap thermal response were 
investigated as part of this study. From the analysis, the following conclusions were drawn: 

a. Convective heat transfer, by itself, may not account for the temperature rise and melting 
during booster cap testing. Further computational and/or empirical quantification of gas 
temperatures and heat transfer coefficients are needed to draw a definitive conclusion as 
to whether convection alone is sufficient. 

b. Zirconia deposition and the subsequent phase change from liquid to solid state may assist 
in booster cap heating and subsequent melting but do not produce booster cap 
temperatures in agreement with the booster cap transient temperature response observed 
during testing. Larger quantities of zirconia deposition increase the propensity to melt 
and accelerate the temperature rise of the booster cap bottom. 

c. Unbumed ZPP deposition can liberate sufficient energy to locally melt through the 
booster cap as indicated by detailed thermal analysis of a hemispherical globule of ZPP 
with a radius of 3.6x10' in and a mass of 5.2x10' mg (1.14x10' lbm). Deposition of as 
little as 20 percent of the unbumed ZPP (~4.6 mg, or I x I O'" lbm) can liberate sufficient 
energy to melt the entire booster cap. Subsequent two-dimensional axisymmetric thermal 
analysis shows that local melt-through can be accomplished with considerably less ZPP. 
From this, it is concluded that « 4.6 mg ZPP burning in contact with the booster cap is 
sufficient to produce the observed response. 

While heat transferred by convection and the phase change of liquid zirconia deposited on the 
booster cover simulator, the unburned ZPP results are a compelling indication that this 
mechanism could be a viable alternative mechanism. There would presumably be a lot of 
variation in the quantity, size and impingement locations of unbumed ZPP that could explain 
inconsistent performance of the NSIs. 
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Appendix F. Phase I Statistical Analysis and Results 

The statistical approach to Phase I was designed to answer the statement of problem: 

Given a limited number of test articles, what differences can be seen 
between the original Al Y-PCA design and the SS V-PCA one in a number of 
responses? Is there a difference in these responses given a dual 
simultaneous firing of a pair ofNSI over a single-NSI firing? 

The Phase I test was designed using some design of experiments (DOE) methods. The statement 
of problem was clearly developed. The size of the test was governed by resource considerations 
rather than a quantitative statement regarding a confidence level on a difference the team desired 
to be able to find, but the subject matter experts had reason to believe this test would be 
sufficient to show clear and useful differences. They were correct in this assessment, as will be 
seen. 

Due to the availability of parts, the runs were not completely randomized. This may have had an 
effect on the data, though it is not quantified or likely even quantifiable. The results, however, 
are clear enough that the qualitative and at least the rough quantitative conclusions and resulting 
recommendations are correct and useful. 

In addition to the direct response data analysis (generally using analysis of variance (ANOVA) 
techniques), some information and data surrounding conduct of the test, such as Trigger voltage, 
were mined for useful observations and anomalies. 
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F.l Analysis and Results 

A summary of Phase I test analysis and results is shown in Table F-l. 

Table F-l. Summary Table of Response Means and 95 Percent Confidence Intervals on the Means 


Response 

SS Mean 

±95 

percent 

A1 Mean 

±95 

percent 

Notes 

Time of Maximum 
Temperature 

776 

120 

959 

120 


Peak Temperature 

2,290 °F 

431 

1,391 °F 

431 

Ignored Run 9 

Temperature at 800 ps 

2,247 °F 

440 

1,162 °F 

440 

Ignored Run 9 

Time of Port A Pressure 
Rise Start 

144 ps 

5 

137 ps 

5 

Ignored Run 9; included duals 

Port A Calculated Peak 





Ignored Run 9; included duals; 

Pressure 

10,729 psi 

289 

7,796 psi 

1,292 

weighted analysis 

Port A Pressure at 800 ps 

4,768 psi 

406 

3,420 psi 

363 


Time of Port B Pressure 
Rise Start 

-295 ps 


-295 ps 


Difference not significant 

Port B Calculated Peak 
Pressure 

4,779 psi 

313 

1,895 psi 

283 

Included Duals 

Port B Time to Peak 
Pressure 

-830 ps 


-1,400 ps 


Runs 9 and 9B may be outliers 

Port B Pressure at 800 ps 

4,588 psi 

300 

1,637 psi 

269 



The Data. This narrative will refer to the Run Number to signify individual trials. The order in 
which its trials were run is different from the Run Number. Run 9 showed some unusual 
characteristics, including being the only A1 run to bum through the diaphragm simulating the 
booster charge, and was rerun as Run 9B. Run 9B ended up as a misfire and was not rerun 
again. It is possible that no ran after Run 3 is valid, but the conclusions presented below do not 
depend heavily on these later runs. 
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Run 

Number 


Order 

Trials Material 
Run 


Number 
of NSIs 


TIME 

OF 

PRESS 

A 

START 


PRESS 

A 

CALC 

PEAK 


TIME TIME 

TIME OF OF PRESS OF 
PRESS PRESS B PRESS 

A PEAK B PEAK B 

START PEAK 


DIA- TIME 
PHRAGM OF 800us 
MAX MAX PRESS A 
TEMP TEMP 


800us 800us 
PRESS TEMP 
B DIA 


2 

1 

SS 

Single 

140 

10926 

156 

255 

4570 

824 

2823 

835 

4638 

4527 

2836 

6 

2 

SS 

Dual 

152 

10419 

168 

151 

9735 

168 



9252 

9196 


7 

3 

SS 

Single 

147 

10781 

166 

324 

4996 

787 

1967 

821 

4933 

4979 

1957 

14 

4 

SS 

Single 

144 

10436 

159 

314 

4441 

870 

1885 

726 

4690 

4232 

1842 

15 

5 

SS 

Single 

147 

11085 

169 

259 

4805 

845 

2485 

721 

4811 

4614 

2352 

4 

6 

Al 

Single 

143 

6881 

157 

286 

1650 

1241 

1308 

1053 

2938 

1580 

1133 

8 

7 

Al 

Single 

138 

8199 

154 

297 

1860 

1480 

1412 

913 

3337 

1307 

1264 

3 

8 

Al 

Single 

130 

8001 

156 

319 

2269 

1440 

1691 

1064 

3857 

1716 

1404 

9 

9 

Al 

Single 

155 

6560 

160 

292 

1848 

1130 

2902 

804 

3077 

1762 

2894 

1 

10 

Al 

Dual 

149 

6456 

164 

149 

6244 

163 

719 

10210 

4863 

4857 


9B 

11 

Al 

Single 

137 

9445 

154 

307 

2149 

1445 

1154 

2875 

3892 

1820 

846 


All PCAs had the same channel area, though the volume was likely different between the A1 and 
SS blocks due to the design differences. 

Dual Firings. Both dual firings resulted in peak temperatures below the temperature required for 
booster charge ignition and times to these maxima. Pressures in these were high. This was a 
clear indication that dual, simultaneous firing of redundant NSIs could result in a chain of events 
resulting in lower reliability instead of the higher reliability intended. 

Parameters Involved in NSI Firing. No important conclusions were drawn from the 
measurements discussed in this section, but the reader should be aware of the observations. 
Traces for Runs 1 and 8 were different from the others. 

Firing Voltage. No important conclusions were drawn from this measurement, but the reader 
should be aware of the observations. An exponentially weighted moving average (EWMA) 
(Roberts, 1959) was used in Figure F-l to smooth the data for visual analysis. 
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Two of the runs, Run 1 (dark grey) and 8 (dark blue), exhibited unusual firing voltage signatures 
over time compared to the other runs. Run 2 was the A1 dual firing. Run 8 was another A1 run 
that seemed to have no other distinguishing responses or features. The apparent differences 
showed up as a noise event at just over 700 ps. Discussion was held on these anomalies, but as 
the firing voltage appeared only to be an issue for the extremely short interval needed for 
actually firing the NSIs and no anomalies were seen during that time frame for these trials, it was 
not explored further. 
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Figure F-l. Phase I Firing Voltages, Smoothed Data, All Runs, Long (6 ms) Time Scale, showing 

Range of Voltages Measured 
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Figure F-2. Firing Voltage, Smoothed Data, All Runs, Time Scale Showing Detail of NSI Firing 

Event, Showing Range of Voltages Measured 


Using a small time scale, visual analysis of this data showed a number of differences between 
trials. Again, the SS runs (run first) appeared to show different features than the Al runs. This 
may have been due to some change in the firing circuit or some difference in the blocks 
themselves, but it was decided that this was only of minor interest. 

Run 9B (dark pink), which held considerable interest or analysis, gave an unusual firing voltage 
trace including comparatively low noise features during the 90-150 ps time interval. 

All the SS runs showed a tight and slightly lower firing voltage than the Al runs, which showed 
some variability. Again, no important differences were seen during the actual time when the 
voltage was required for NSI firing, so this was not extensively explored. Again, this 
information is provided in the interest of completeness. 

Firing Current. It is in this trace that Runs 1 and 8 appear to have experienced shorts or some 
similar disturbance. Again, Run 1 is the dark grey trace and Run 8 is the medium blue one. 

Both show their differences after the NSIs actually fire. The event at 5,000 ps is discontinuance 
of energy to the NSI. 
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Figure F-3. Firing Current, All Runs, Long (5 ms) Time Scale, Showing Range of Amperages 

Measured 
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Figure F-4. Firing Current, All Runs, Time Scale Including NSI Initiation through Firing Event, 

Showing Range of Amperages Measured 
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Temperature. The IR setup, again, did not register temperatures below 572 °F. 

Temperature traces indicated differences between runs. Most of the SS trials showed a higher 
peak temperature than the A1 trials. 
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Figure F-5. Temperature Readings, Single Firings Only, Over Time Scale Including NSI Initiation, 
NSI Firing, Temperature Registration on IR Instrument and an Approximate Millisecond after 

Temperature Rise 

Run 9B’s initiator clearly fired, but the temperature climbed only slowly and never reached a 
high value. This could have resulted in a misfire in actual use. Remember that this was not a 
dual firing. It is possible that a redundant NSI would have resulted in a successful booster firing. 
It is also possible that the firing circuit was to blame and a second NSI would have experienced a 
failure due to a common cause. 

If it is assumed that the possible misfire was due to a discrete event that could occur in actual use 
conditions, it cannot be assumed conclusively that the event cannot happen in SS PCAs. There 
are too few trials to be able to conclude this. However, if one assumes that the 9B run is part of a 
continuum of responses that Al Y-PCAs could give, which is reasonable given this graph 
(Run 9B looks qualitatively like a low-response cousin to Runs 3, 4 and 8), it seems reasonable 
to believe that the SS design will be less likely to produce this type of outcome. 

Peak Temperature. If Run 9, the only Al run having diaphragm burn-through, is included, the 
team cannot show, using the data, that there is a difference between Al and SS designs. 

Removing Run 9 allows a conclusion that SS results in a higher peak temperature than Al at 
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greater than 95 percent confidence. (From this point on, it will be assumed that confidence is 
greater than 95 percent unless otherwise indicated and the confidence level will not be generally 
stated explicitly.) Run 9B was left in this dataset for analysis. Its removal results in little 
difference in conclusions. 
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Figure F-6. Diaphragm Peak Temperature, Data Including Run 9 
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Figure F-7. Diaphragm Peak Temperature, Data Not Including Run 9 

Time to Peak Temperature. Quantitative analysis shows clearly that time of maximum 
temperature is lower for SS than for A1 PC As. 
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Figure F-8. Time of Maximum Temperature from Time NSI fired, Run 9B Not Included 
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Figure F-9. Time of Maximum Temperature from Time NSI Fired, Runs 9 and 9B Not Included 
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w/ Runs 9, 9B 

Figure F-10. Time of Maximum Temperature from Time NSI Fired, Runs 9 and 9B Included 

If Run 9B is not included in analysis, it can be concluded with greater than 95 percent confidence 
that SS gives shorter Time than Al. Adding 9B back in will increase the difference, as it has a 
long Time to Max Temperature. However, because 9B increases spread tremendously, it has the 
paradoxical effect of decreasing confidence in the hypothesis. Removing both 9 and 9B 
increases confidence in this statement. Note that the analysis without Runs 9 and 9B makes 
Run 4, the first Al trial run, look like a potential wild point. 

Diaphragm Temperature at 800 us. Again, Run 9 affected results; including Run 9 resulted in a 
finding of no significant difference, while leaving it out resulted in rejecting the hypothesis that 
there was no difference between SS and Al PCA designs. 
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Figure F-ll. Diaphragm Temperature at 800 ps, Run 9 Included 
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Figure F-12. Diaphragm Temperature at 800 ps, Not Including Run 9 

Pressure. Port A Measured Pressure. In visual analysis of the traces, pressure reached a peak 
close in time for both designs. A small blip occurred just before pressure rise to the peak 
pressure. Stainless PCAs resulted in higher pressures after peak than for the Al design. Ringing 
oscillations were seen and showed a reasonably distinct signature. Oscillations damped out by 
800 ps. 
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Figure F-13. Traces of Port A Pressures 


Time of Pressure Rise Start, Port A. There may have been a difference between SS and A1 
designs, but it was small and borderline statistically significant. Not including Run 9 and 
including dual/single firing as a factor improved resolution and resulted in a statistically 
significant difference seen between A1 and SS. It is not obvious that the 6.6 ps difference seen 
between A1 and SS is significant in engineering terms. This parameter does not have much 
practical significance and is probably influenced by minor differences in how the grease is 
packed into the pressure transducer sense ports. 
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Figure F-14. Time of Rise in Pressure, All Runs 
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Figure F-15. Analysis Showing Difference between Dual and Single Firings and SS and Al Firings If 

Particular Data Are Included (See Text) 


Port A Extrapolated Peak Pressure. A team member (W. Sipes) developed a method to calculate 
an approximate peak pressure for each test that eliminated response issues in the data. The 
method is described elsewhere in this report. This data showed a significant difference between 
SS and Al (p = 0.003, or confidence at 99.7 percent). Graphical results are shown in the 
Figure F-16. 
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Figure F-l 6. Peak Pressure in A 


However, there was a potential issue with variability growth toward end of testing. The 
difference in variability between A1 and SS is shown in the following graphs. Figure F-l 7 shows 
the difference in variability. 
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Figure F-l 7. Residuals (Actual Data Value — Predicted Value) Compared between Al Y-PCA 
(Left-Hand Points) and SS V-PCA (Right) for Peak Pressure in A 
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Figure F-18 shows the increase in variability as testing continues. Cook’s distance 
http://data.princeton.edu/wws509/notes/c2s9.html is a measure of how far a (standardized) 
datum is away from expectation given a model. 
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Figure F-18. Analysis of Residuals; Bottom Axis Is Order of Trials Run; Suggestion Is that Variability 
Increased as Testing Progressed; Issue Would Have Been Able To Be Separated from Other 
Conclusions Had Order of Tests Run Been Randomized 
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Figure F-19. Peak Pressure in A: Comparison Using Weighted Least Squares 
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Port A Pressure at 800 jus. The difference in variability appears also to occur at long times. The 
following analysis uses simple regression, as the difference in variabilities between A1 and SS is 
lower at this time scale. If weighted least squares had been used, the difference in means (and 
the means estimates) would have been the same, but the confidences on the means would have 
been different: larger for Al, smaller for SS. 
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Figure F-20. Pressure in Port a at 800 /. is 


Port B Measured Pressure. Measured pressures in single NSI firings for the SS design were 
easily seen to be higher than for Al PCAs. The traces for each were highly characteristic and 
clearly showed that the SS design allowed much higher flow rates across the top of the booster 
charge and into the opposite NSI cavity. 
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Figure F-21. Pressure Traces for Non-Fired Port 


Time of Start of Pressure Rise , Port B. The difference between SS and Al for the time at which 
pressure began to rise in the non-fired port was not significantly different statistically. Note that 
this does not mean that there is no difference, only that the difference could not be discerned 
with this data given its size and variability. 
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Figure F-22. Time of Start of Pressure Rise 


Port B Pressure at 800 jus. Again, as is abundantly clear in the pressure traces, Al pressure in 
the non- fired port is considerably lower than in the SS non- fired port at late time slices. 
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Figure F-23. Port B Pressure at 800 jus 


Difference between Port A and Port B Pressure. It takes some time for pressure to equilibrate 
between the fired and non- fired ports. A couple of features are interesting in the following trace 
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of this value. For SS NSIs, the difference is larger than for A1 during the time from just after the 
initial peak to about 500 ps. It then decreases quickly, and for a short time, the Pressure in A 
appears to be lower than that in the non-fired port for a few tens of a microsecond. A1 PCAs 
may ring for a bit longer than the SS and the ports take more time to equilibrate than for SS. It is 
also possible that there are two families of A1 firings in this data: Runs 9 and 4 have similar 
traces, while the others overlap each other at a parallel, higher value as compared to Runs 9 and 
4. 
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Figure F-24. Trace Showing Difference in Pressure between Ports A and B Over Time 
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Appendix G. Phase IIB Statistical Analysis and Results 

Phase IIB was designed to evaluate the relationship between firing skew and the cross-sectional 
area of the flow passages and various responses relating to factors involved with firing a booster 
charge. It was assumed that the following model would apply. 

Y = Po + Pi Skew + P 2 Area + P 12 Skew x Area + Pn Skew 2 + P 22 Area 2 

where the P’s are linear regression parameters fit using the data. 

G.l Analysis and Results 

The analysis results were clear: within the tested range of the two inputs studied, no-fires due to 
insufficient temperature at the booster interface are entirely possible at low skews, even in 
designs having enlarged channels. A low-temperature response occurred in this test at conditions 
of high skew and channel area (Run 3). Because of this, it cannot be ruled out that a no-fire at 
more than 500 ps could occur, but the team believes it much more likely that the response is due 
to installation of two booster simulator diaphragms in error. In fact, the inner of the two 
diaphragms was burnt through in that trial. Temperatures increase with increasing skew but are 
fairly insensitive to channel area. Channel area affects pressure and increases variability of the 
time it takes to attain peak temperature. 

The data analyzed is shown in Figure G-l . It was developed from the raw traces and other 
sources by team members. 
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Figure G-l. Data 


Channel Area was used as the input factor for these analyses rather than Channel Diameter. 
Actual Skew was used in place of the intended value. 
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The models created below are not primarily physics-based and are not intended to be. The exact 
shapes of the response surfaces created, particularly in low-skew conditions, cannot be, nor are 
intended to be, strictly representative of the physics. They provide a useful predictive picture of 
what is happening. If greater fidelity is required, the analyses and response surfaces point the 
way to where follow-on tests will be most informative. 

Temperature. Two clear misfires, Runs 1 and 3, were observed: Run 1 did not result in 
temperatures above the infrared (IR) measurement threshold of 574 °F and Run 3 barely 
exceeded the threshold and took a considerable amount of time to do so. These were excluded 
from analysis in most cases. Run 15, another apparent misfire, was also often excluded, as 
shown in the figures below. It may be more appropriate to perfonn an analysis where these 
points are censored. Details will be discussed in each case below. 

Overall Maximum Temperature was influenced by Runs 1 and 3. These misfires clearly 
include a different heat transfer mode than most of the other runs. These were removed from the 
data for analysis. A censored analysis was not performed: because of the way the data were 
censored, there would be no effect on the analysis conclusions. 

Run 15 also showed a low Peak Max Temperature of 699 °F. This run also is probably 
dominated by an undesirable heat transfer mode. Analysis including this data point suggests it 
may be since it shows as a possible wild point using Cook’s D test, but not any other diagnostic 
measures used. Analysis with Run 15 included results in the following model: 

Overall Maximum Temperature = 

+ 1485. 

+3.359 * Actual Skew 

Given Run 15 is not unusual since Overall Max Temperature increases predictably with 
increasing Skew, Channel Area does not have a discemable effect. 
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Overall Max Temp 
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Figure G-2. Maximum Temperature, Prediction Response Surface Given Various Actual Skew and 
Channel Area, Data Including Run 15 (Hollow Point Marker) 


Analysis without Run 15 shows a weak (p = 0.08) and slight effect of Channel Area. It also 
suggests an increasing rate of increase in Max Temperature with Skew. 


Overall Maximum Temperature = 

+ 1629. 

+42670 * Channel Area 

-0.6889 * Actual Skew 

+0.007156 * Actual Skew 2 
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Design- Expert® Software 
Factor Coding: Actual 
Overall Max Temp 
* Design points above predated 
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Figure G-3. Maximum Temperature, Response Surface Created Ignoring Run 15 


Analysis of Diaphragm Temperature at 1,600 ps did not include Runs 1 and 3. An analysis using 
censored data would not change the analysis results. Run 15 did not flag as an outlier in tests of 
residuals. 


Diaphragm Temperature at 1,600 ps = 
+936.8 

+3.452 * Actual Skew 


NESC Request No.: 09-00596 


© 

NASA Engineering and Safety Center 
Technical Assessment Report 

Document #: 

NESC-RP- 

09-00596 

Version: 

1.0 

Title: 

Pyrovalve Booster Interface Temperature Measurement 

Page #: 

61 of 92 


Design- Expert® Software 
Factor Coding: Actual 



Figure G-4. Diaphragm Temperature at 1,600 ps 


Run 14 flags a possible unusual point using a Cook’s D statistic, but it has substantive effect on 
predictions and none on conclusions. It was a replicate run of Run 3, but it had the correct 
number of diaphragms installed (1). 

Time of Maximum Temperature was difficult to interpret. Run 1 resulted in no data as its 
temperature never ventured above the detection threshold of the IR unit. Values for Run 3 and 
Run 15, the other two runs that never achieved 1,000 °F, were also quite high. 

The data could not be modeled as-is using ordinary least-squares regression assuming normally 
distributed residuals, even if Runs 1 and 3 are left off. However, if Run 15 is also ignored — the 
assumption then being that only the trials reaching at least 1 , 1 00 °F are modeled — a useful 
relationship can be modeled. 

A log transformation of the response was applied to account for an increasing variability with 
increasing response value. It also has the beneficial effect of making negative values impossible, 
which is useful for constructing predictions and uncertainty bounds on them. This model fit the 
data reasonably well; Run 2 was somewhat high compared to expectation, but it did not flag as a 
clear outlier. The model: 
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Logio(Time of Max Temp) = 

+2.340 

+4.492E-003 * Actual Skew 

-6.548E-006 * Actual Skew 2 


The graph of this response surface is in Figure G-5: 


Design-Expert® Software 
Factor Coding: Actual 
Original Scale 
Time of Max Temp 

• Design points above predicted value 



Figure G-5. Time of Maximum Temperature, Model Formed Using Data for which Maximum 
Temperatures Reached at Least 1,100 °F. Run 2 Is Signified by a Hollow-Center Marker. 


Removing Run 2 provides a better fit all around, but the model and predictions change little. It is 
apparent from the data that at low skew, maximum temperature is reached quickly compared to 
higher-skew trials, or near 200 ps at 0 skew. It increases to a median value just under 
1,000 ps at or before 235 ps skew and appears to plateau. Variability appears to increase with 
increasing skew. Again, this assumes peak temperature achieves 1,100 °F. Channel area was 
not shown to have significant effect on this response. 

If only the data for trials run at medium and high factor settings are examined, excluding Run 3 
for skew, a significant linear regression model could not be fit. That is, neither skew nor area 
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had a discemable effect on Time to Peak Temperature when Skew was above about 235 us. The 
data can thus be described as randomly distributed. A reasonable fit to a nonnal distribution was 
made when logio ( Time to Peak Temperature) was modeled. The expected mean was 1,040 ps, 
with 95 percent confidence limits on that mean estimate of 800, 1360. This model will be used 
in predictions. 

Time to 1,000 °F excludes Runs 1, 3 and 15 because the temperature never reached that level. A 
model suggested by the data resulted in the following response surface, shown in Figure G-6. 


Design-Expert® Software 
Factor Coding: Actual 
Time to 1000F 

• Design points above predicted value 

• Design points below predicte 

XI = A: Channel Area 
X2 = B: Actual Skew 



channel Are; 


Figure G-6. Time to 1,000 °F 

The surface was not satisfactory to the assessment team. It is possible that there are two 
processes: one that drives a quick, 200 ps jump to 1,000 °F; another that results in a more 
leisurely pace of 500-600 ps. It is also quite possible that above 0 skew, time to 1,000 °F is 
distributed randomly 

Calculated Peak Pressure in the A and B Channels. This was the highest pressure recorded in 
the A (first-fired) and B (second-fired) channel. Simultaneous firings result in no first-second 
ordering, but two pressure readings were taken. Some estimation was used due to noise and 
overshoot in the data. 
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This data was relatively straightforward to analyze and fit two linear response surfaces. Run 14 
showed as a possible wild point in the Peak B response. There is evidence that this run’s 
window leaked slightly, which jibes with the Peak B analysis results. Removing the point does 
not change the overall conclusions. The point was left in the model. 

The prediction model fit for Peak A is shown below. 


Press A Calc Peak 
+10380 
-102700 
+9.152 
-2933 
-183500 
+148100 


* Channel Area 

* Actual Skew 

* Channel Area * Actual Skew 

* Channel Area" 

* Channel Area 2 * Actual Skew 


Fit to assumptions was good although the goodness-of-fit of the model was only fair. Removing 
Run 14 did not change this and extensive diagnostics and model-fitting were not performed to 
find the source of the issue. Again, the overall conclusions due to this data are clear, and it is 
expected that a better model fit would only result in more accurate point predictions without real 
value to the aims of this test. 

A 3 -dimensional view of the fit including the data points is shown in Figure G-7. Run 14 is 
highlighted by an open marker. 
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Figure G-7. Peak Pressure in First-Fired Chamber / Port A 


Highest Peak A Pressure occurred at the most-commonly designed low Channel Area and higher 
Skew times. Note that Channel Area produced only a —1 ,000-lb drop from the minimum to 
maximum tested Area on Peak A Pressure at low skew. Similarly, at low channel areas, Peak A 
Pressure increased only about 1,000 lbs from 0-500 ps skew. Both look linear at these settings. 
However, with large Channel Area and Skew, Peak A Pressure falls off noticeably. There is a 
clear curvature in the response to Channel Area. 

Peak B data was fit to a model that was in some ways similar. 


Press A Calc Peak 
+10510 
-83450 
+21.60 
-649.6 
-0.01976 


* Channel Area 

* Actual Skew 

* Channel Area * Actual Skew 

* Actual Skew 2 
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The response surface described is shown in Figure G-8. 
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Figure G-8. Peak Pressure in Second-Fired Chamber / Port B 

Again, at no skew, Channel Area had little effect on Peak B Pressure. Peak B Pressure rises 
from here with higher skew. The model appears to show a peak around 300 ps skew, but it 
seems possible that given the physics, Peak B Pressure levels off at a pressure that is a function, 
which decreases with increasing Channel Area. In either case, the Peak B Pressure peak value is 
highest at low Channel Area. 

Figure G-9 shows the relationship between the two response surface models. The lower surface 
is Peak A Pressure and the upper is Peak B Pressure. 
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Desagn-Expertt&SoftraB 

Factor Codtog: Actual 
Press B Calc Peak 



Figure G-9. Overlain Plots of Peak A (blue / lower) and B (black / upper) Pressures. Point Markers 

with Hollow Centers Are Run 14 


As expected, the modeled surfaces meet at the 0-skew point. The Peak Pressure after B fires 
builds predictably above the Peak A, again an expected result. Larger channel area can only 
decrease Peak Pressure given any Skew. Given the current default 0.060-inch-diameter channel 
design, Peak Pressure after the first NSI fires is stable and predictable, and pressure builds 
considerably onffs firing. 

Pressures at Key Times. Analysis of Pressure in Channel A at Time of Channel B Pressure Peak 
data results in a fairly complex surface. It was modeled as follows. 


NESC Request No.: 09-00596 


© 

NASA Engineering and Safety Center 
Technical Assessment Report 

Document #: 

NESC-RP- 

09-00596 

Version: 

1.0 

Title: 

Pyrovalve Booster Interface Temperature Measurement 

Page #: 

68 of 92 


A Press at B Press Peak 


+10780 

-225900 * Channel Area 

+2.332 * Actual Skew 


-5207. * Channel Area * Actual Skew 


+9345000 * 

-0.01024 * 
+151700 * 

+5.470 * 


Channel Area A 2 
Actual Skew A 2 

Channel Area A 2 * Actual Skew 
Channel Area * Actual Skew A 2 


This model provides a close fit to the data, which show little variability at each factor setting 
combination. It appears that the linear drop of -1,000 lbs is seen again from low to high 
Channel Area at 0 Skew. Increasing skews, as might be expected, produce a steeper, but still 
linear drop from -10,000 lbs to -6,000 lbs over the range 0-5,000 ps at lowest area. Again, this 
pressure drops markedly from these cases. It can be conjectured that the dip in the response 
surface just before the high-Skew/high-Channel Area factor combination is an artifact of the 
modeling process, and the true process may in fact level off at approximately 4,000 lbs in this 
region. 
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Figure G-10. Pressure in First Chamber at Time of Peak Pressure in Redundant NSI 9 s Chamber 
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The Difference in Pressure at the Time of Peak Pressure in B shows little if any dependence on 
Channel Area and a markedly decreasing dependence on Skew. A floor at approximately 
-1 1,000 lbs seems reasonable. The data suggest a slight downward bow with minimum due to 
Channel Area with minimum at ~ 0.0070 in 2 , but it is small. This tenn was arbitrarily left out of 
the model. 


A Pressure at Time of B Pressure Peak 
-90.96 

-44.02 * Actual Skew 

+0.05412 * Actual Skew A 2 
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Figure G-ll. Difference between Port A and Port B Pressures at Time of B Peak Pressure 

Pressure in Channel A at 1,600 jus depends on Channel Area but not appreciably on Skew. 
Clearly the effect of Skew has diminished by this time and total PC A volume and/or metal area 
available for heat dissipation drive this. 
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Pressure in Channel A at 1 600 ns 
+9942. 

-430000 * Channel Area 

+ 19150000 * Channel Area 2 


At the default Channel Area, one could expect pressures of around 8,800 lbs. It may decrease to 
a floor of around 7,500 lbs. 
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Figure G-12. Pressure in Port A at 1,600 jus. Run 14 Is the Hollow-Center Point in the Graph. 


Time of Peak Pressure A decreased slightly but fairly dependably with increasing Channel Area. 
Run 6A attained peak pressure at an unusually low time, according to this model, but the test 
value (212 ps, blue-circled point in the graph) was not far off the expected value (221 ps) in real 
terms due to a fairly small spread. Nevertheless, it does look unusual in the plot of the response 
surface. 
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Figure G-13. Time of Peak Pressure in Port A 


The model fit is as follows. 

Time of Press A Peak = 

+224.3 

-528.3 * Channel Area 

No linear model was found to fit Pressure before B Begins to Rise. This measure is not defined 
for 0-skew runs. The natural logarithm of the data for the remaining domain can be fit to a 
normal distribution. This would suggest that the variability in the Pressure before B Rise data 
has little to do with either Area or Skew and is due to random (uncontrolled) factors, and its 
variability can be modeled using a lognormal distribution. 

The QQ plot shown in Figure G-14 indicates the fit to the distribution. A perfect fit of the data 
to the distribution would have the points lying exactly along the centerline. Note that it takes a 
large dataset to pin down a distribution, so this information is presented as being useful, but not 
certain. 
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Quantile-Quantile Plot 



Figure G-14. Pressure before Port B ’s Pressure Begins to Rise, Modeled Using a Lognormal 

Distribution 


The distribution shown in Figure G-15 is the probability density function (PDF). 
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Density Function 



Press before B Rise 


Figure G-15. Probability Density Function for Distribution Best-Fit to Data and Indicated in 

Figure G-14 


Its parameters: 


Lognormal 


mean = 3683 


standard deviation =1017 


Log scale: mean= 8.175 


Log scale: std. dev. = 0.271 1 
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Time of Pressure B Peak is clearly due to Skew, which because skew is directly related to B's 
firing time, was not a surprise. A case can be made that there is lower variability at lowest 
Channel Area than at other areas tested. Skew had no effect. 

Size of Hole Made in Diaphragm 

Responses on a percent scale are bounded at 0 percent and 100 percent. These are best handled 
using a transformation of the response that reflects this such as a Logit but, in practice, if 
responses are in the center of this range, this may not be necessary. After confirming that there 
is no substantial difference in conclusions reached by the more complex model and a model 
reached using no transformation, the untransformed response was used in analysis. 

This measure’s analytical results depend on how no-fires are handled. No-fires reliably resulted 
in 0-percent diaphragm hole. All of the 0-skew runs resulted in 0-percent puncture area. The 
one no-fire at the largest skew setting, Run 3, presents a problem, particularly because the team 
was not able to collect direct data from its replicate point, Run 14. There is a sufficient still from 
Run 14 that clearly shows substantial diaphragm burn- through. 

There is no reason to call Run 14’s response 0 percent; in fact, it may have the largest burn- 
through of any test. It could be treated as a point censored at 0 percent, the most correct method, 
but doing so would result in Run 3 driving the analysis. Censoring at a higher value would work, 
but was not done at the time of this report’s writing. 

It seems clear that Run 3 is somehow anomalous. One method is to ignore or right-censor this 
point. Censoring this value has little effect on results and presents a challenge, which takes time 
to solve, so it was decided to ignore the point. A placeholder value of 50 percent was used for 
Run to avoid the previously mentioned analysis time issue. Using 70 percent results in an 
anomalously high value, an outlier; it also suggests an effect due to Area. 

It is also possible that these assumptions are entirely incorrect and that the large S&ew/large Area 
point simply results in unpredictable burn-through. 

A fit model, assuming Run 3 does not exist and Run 14 is 50 percent: 

Diaphragm Fract. Hole, Nonzero Skew = 

+0.1427 

+0.0005861 * Actual Skew 

The response surface (Run 14 is the hollow marker) is shown in Figure G-16. 
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Design-Expert® Software 
Factor Coding: Actual 



Figure G-l 6. Size of Hole in Booster Simulator Diaphragm 


Clearly, low Skew could result in 0 bum- through, an undesirable result. Given Run 14 is due to a 
special cause or given a no-fire does not occur, it appears that larger Skew results in larger 
fraction of the diaphragm burning. 

Phase IIB Confirmation Run . A confirmation run is a test trial withheld from the dataset, which 
had been used to construct predictive models. The models are created without this ran, but it is 
then used to check the reasonableness of the model predictions by comparing the results of the 
withheld test to the predictive models. One useful way to evaluate this comparison is through 
the use of prediction intervals (Montgomery, Peck and Vining, 2001: Introduction to Linear 
Regression Analysis, 3 ld Ed., Wiley Interscience, New York, NY). Prediction intervals are test 
result ranges within which one would expect a future observation to lie with some specified 
probability, given a model fonned from data that has already been observed. This is, of course, 
exactly what is desired here. 

A single confirmation ran at factor settings Nominal Skew = 500 ps, NSI flow channel cross- 
sectional area = 4 times the nominal (Channel Diameter = 0. 120 inches) was performed. This 
point was chosen because of the issues with Run 3, inadvertently ran with two booster simulator 
diaphragms instead of one, and Run 14, having indications of a pressure leak. The results were 
as follows: 
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Temp 

Run 

3A 

0.120 
4 times 
nom) 

488 

8243 

214 

>3630 

(above 

limit 

ofIR) 

299 

4140 

10,114 

700 

7773 

2958 

All times are measured from start of firing signal. 


Note that the maximum temperature reported is at the upper limit of the pyrometer. The 
pyrometer sensor was saturated for some period of time, so the actual maximum temperature is 
higher than the value reported, but it cannot be used. 

To evaluate reasonableness, 95-percent prediction intervals (PI) on the models created as earlier 
described were compared to the values above. The summary is shown in Table G-l . 

Table G-l. Confirmation Run 3 A Compared to Model Predictions and Prediction Intervals 
Confirmation Report 

Two-sided Confidence = 95% 


Response 

Name 

Level 

Low Level 

High Level 


Run 

3A 

Channel Area 
Actual Skew 

0.0123 

0.0028 

0.0123 


0.0123 

488 

0 

486 

488 US 

Prediction 

Std Dev 

SE (n=l) 

95% PI lower 

95% PI upper 

Run 

3A 

Overall Max Temp 

3522 

269.0 

331.3 

2783.7 

4260.0 

>3,630 

1600us Temp Dia 

2622 

286.1 

313.0 

1945.3 

3297.8 

2,958 

Time of Max Temp 

1040 



416 

2607 

299 

Press A Calc Peak 

6876 

380.1 

459 5 

5864.9 

7887.8 

8,243 

Time of Press A Peak 

217.8 

3.39 

3.68 

209.9 

225.7 

214 

Press Before B Rise 

not analyzed 

4,140 

Press B Calc Peak 

11417 

712.8 

863.5 

9535.4 

13298.2 

10,114 

Time of Press B Peak 

not analyzed 

700 

1600us Press A 

7550 

438.8 

480.7 

6519.4 

8581.4 

7,773 


For most measures where a model was constructed, Run 3A falls within the 95-percent PI upper 
and lower bounds. This gives confirmation that these models are useful for prediction of future 
test runs. Note that prediction intervals are not intended to predict the range within all future test 
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runs can be expected to lie. That construct is known as a tolerance interval and is understood 
and calculated quite differently. Tolerance intervals are normally considerably wider than 
prediction intervals. Maximum temperature was not evaluated, as its true value exceeded the 
upper limit of the IR unit’s capability. 

For the case of Peak Pressure in Channel A and Time of Maximum Temperature, Run 3A’s 
measured value is outside the 95-percent PI. By chance, one might expect one out of every 20 
observations (1-0.95, or 5 percent) to be outside these limits. However, given six (not including 
Max Temperature ) independent measurements to be compared at one time, there will be a greater 
than 5 -percent chance that any one of these measurements will be exceeded. The chance that 
there will be at least one exceedance in six of these intervals simply due to chance can be 
calculated using the Bonferroni correction (Montgomery, Design and Analysis of Experiments, 
6 th ed., Wiley, 2004), which comes out to more than 26 percent. Clearly, it is quite possible to 
experience at least one exceedance of these limits in six chances to do so without needing to 
resort to any explanation regarding special causes. 

To reduce the chance that a single exceedance of these six prediction limits would be observed 
due to Run 3A to 5 percent, the confidence level to 99.17 percent would need to change. Doing 
so increases the upper limit on the prediction limit (now a two-sided upper 99.17-percent 
prediction limit) to 8713, bounding 3A’s value. The Time of Max Temperature two-sided lower 
99.17-percent prediction limit is 271 ps, again bounding 3A’s observation; 3 A appears to be 
bounded by the models created, meaning that this confirmation run does not invalidate these 
models. 

It is possible that the model predicting Peak Pressure in A does not well represent high 
skew/large area. Run 3A’s Peak Pressure in A compared to the predicted model is pictured in 
Figure G-17. The point included in that corner is Run 3; the model was constructed without 
Run 14. It is possible that the response surface model is not adequate in that corner of the input 
domain and Run 3A is reasonable. This would imply that Peak Pressure levels off more sharply 
at long Skew and large Channel Area than modeled, perhaps near 8,000 psi. 
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Design-Expert® Software 
Factor Coding: Actual 
Press A Calc Peak 



Figure G-l 7. Response Surface Model, Peak Pressure in Channel A. The Green “<8>” Marks Run 3 A ’s 
Value (8243 psi). Run 3A Was Not Considered when this Model Was Built. 


The confirmation run appears to underline the robustness of this test planning using DOE, the 
test performance, and the models created. It does not appear to indicate that most of the models 
inadequately predict their respective responses. It does not appear to call into question the 
findings, observations and NESC recommendations of this team. 

Analysis Excluding Runs 3 and 15. An analysis of the data was performed without either of the 
data points collected at 500 ps and 0.0123 in“. It was expected that, since the test matrix was 
constructed using DOE methods, and most of the data and analyses were reasonably 
straightforward and well-behaved, this would have minimal effect on conclusions. 

In fact, this was the case. The temperature cases’ response surfaces were barely affected by 
removal of Run 15; Run 3 had already been removed before analysis. Run 14 was no longer 
flagged as a possible outlier in the 1 ,600 us temperature analysis. Time of Maximum 
Temperature was more clearly more variable at higher Channel Area. 

Peak pressure predictions were affected somewhat. The surfaces, now not tempered by the high 
Skew/ high Channel Area points, were slightly less complex. Given the values measured at this 
high / high point, however, it seemed clear that the surface was extrapolating in that corner of the 
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response surface, whereas predictions at that corner still matched the actual, but missing, 
temperature measurements. The models and surfaces describing Pressure at 1,600 ps and Time 
of Peak Pressure in A looked similar whether or not the high/high data was included. The size of 
the hole in the diaphragm was also little affected by presence or absence of these points. 

Optimization and Summary 


A summary of effects of factors and runs that appear to give the individual model problems is 
shown in Figure G-18. 


CL 

E 

0 ) 


Response 


Overall Maximum Temperature 
Diaphragm Temperature at 1600 ps 

Time of Maximum Temperature 

Time to 1000F 


Peak Pressure in A (Calculated) 
Peak Pressure in B (Calculated) 
Pressure in A at B Peak Pressure 
A Pressure at Time of Peak B 
Pressure in A at 1600 ps 
Time of Peak Pressure in A 
Pressure before B Pressure Rise 
Time of Peak Pressure in B 


Hole Size of Hole Made in Diaphragm 


Significant Effect of 
Factor 



Issues 


Run 1 


Run 15 


Flat response skew > 230 ps 
Runs 1 , 3, 15 (Max T< HOOF) 


No model 


Run 14 


Run 6A 


No model 


Run 3; Run 14 ; 0 skew 


Figure G-18. Summary of Significant Effects. A Green Box Indicates a Statistically Significant Effect 

of the Factor on the Response. 


An optimization routine was run. This optimization analysis was not meant to determine 
absolute optimum values for input factor settings in this case. Instead, it was meant to roughly 
determine what input factor settings would produce desirable results and look for issues. 

Team subject matter experts were polled for their estimation of importance of the responses to 
whether the booster would fire or not. There was general, though not unanimous, agreement. 
Individual responses will not be supplied, but Table G-2 shows a summary of the combined 
responses. 
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Table G-2. Subject Matter Expert Responses to Booster Fire 

Lower Upper 

Name Goal Limit Limit Importance 

1600 ps Temp @ Diaphragm 


Overall Max Temp 


Diaphragm Fract. Hole, Nonzero Skew 


A Press at B Press Peak 


Diff in Press at B Peak 


Time of Press A Peak 


A possible graph of the desirability of the total of responses weighted by team opinion over the 
range of the input factors Skew and Area is shown in Figure G-19. 

Design-Expert® Software 
Factor Coding: Actual 
Desirability 

XI = A: Channel Area 
X2 = B: Actual Skew 
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Figure G-l 9. Desirability Function Suggested by Expert Ratings of Relative Importance of Responses 
Combined with Response Surface Models in this Analysis 
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The region was cut off at 100 ps skew to reflect the clear issue of no-fires. This cutoff point is 
arbitrary, and characterization of the probability of no-fires given a skew value is left to future 
testing. The most desirable region of this space is with higher skew. Temperatures under 
1,000 °F were weighted as unacceptable, but greater temperatures were given a high and only 
gradually increasing weight, resulting in a flat desirability over much of the region. 

The analysis of Phase II data clearly shows: 

• Low skew times can result in no-fires. 

• Maximum temperatures are clearly affected by skew, with longer skews associated with 
higher temperatures. 

• The time at which these temperatures are reached is lengthened with larger channels. 

• Designs using larger channel areas should be reconsidered, particularly if the sole aim of 
their design is to prevent no-fires. 

In addition, a no-fire at 500 ps skew indicates a possibility that skews greater than 500 us may 
be necessary and/or that channel areas should not be increased much beyond 0.0028 inches. 
However, it seems most likely that this event actually was due to an incorrectly assembled test 
article. 
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Appendix H. Numerical Simulations of Single and Simultaneous Dual Firing 

Initiators in the SS V-PCA Design 

The objective of the numerical simulation (computer modeling) is to provide a detailed fluid 
dynamic analysis of the single and dual, simultaneous NSI firings in the SS V-PCA design. This 
work did not include the A1 Y-PCA design. 

Some simplifications have been made to allow for a quicker and more efficient preliminary 
investigation. Specifically, the heat transfer to the solid walls was neglected and an inviscid gas 
was assumed. All combustion of the ZPP was assumed to take place in the NSI. That is, 
continuing combustion of any ejected large particles was neglected. In addition, the calculations 
reported here were conducted on a 2-D planar configuration. For the dual, simultaneous firing 
simulation, the difference in firing times of the NSIs (skew) was chosen to be 0 ps. 

H. 1 CRUNCH CFD® Numerical Framework 

The numerical simulations were performed using the CRUNCH CFD R multi-physics numerical 
code [2]-4] that has been developed by Craftech® and is commercially-licensed. This code is a 
three-dimensional, unstructured finite -volume solver for viscous, generalized fluid simulations 
that allow hybrid element grids (i.e., tetrahedral, prismatic, pyramid, and hexahedral cells). 

Some of its key capabilities are: (i) accurate transient shock wave propagation in high-pressure 
gas, (ii) strongly coupled conjugate heat transfer capability to model temperature rise in solid 
material, and (iii) multi-phase capability to model solid propellant particulates as well as molten 
metal that has eroded from the PCA casing. 

H.2 Numerical Results for 2-D Planar Configuration 

For this work, inviscid flow was assumed and combustion of the propellant charges was not 
modeled directly. Instead, a zone of high pressure and temperature was initialized in the region 
of the charges (top of the NSI cavities) as shown in Figure H-l(a), Figure H-2(a), Figure H-3(a), 
and Figure H-4(a). The pressure and temperature levels (82.7 MPa and 2,850 K) were derived 
from an analytical solution of a closed-bomb problem for the complete burn of a typical NSI 
charge. The flow-field is initially quiescent and atmospheric conditions are assumed in the areas 
of the PCA outside of the high-pressure zones. In addition, a trace amount of 20 pm diameter 
ZrC >2 particles were uniformly distributed in the high pressure/temperature zones to track the 
particulate location since heat transfer from the particulate phase is expected to be a major factor 
in the gas burning through the booster cap. 

The pressure contours for the single and simultaneous dual firings at 50 ps, 100 ps, and 150 ps 
area are shown in Figure H-l and Figure H-2. It is clear from these figures that the pressure rise 
in the booster cap region for the dual, simultaneous firing case reaches and sustains a much 
higher peak level than the single firing case due to interaction between the flows from each flame 
channel. In contrast, the pressure in the single firing case gets relief from the unused flow 
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channel, which is at a low ambient pressure. However, the high pressure in the booster cap 
region for the dual firing case effectively creates a stagnation zone and reverses the flow back 
into the flow channels. This effect is more evident when looking at the Mach number contours 
in Figure H-5 and Figure H-6. The single firing case shows a high Mach number in the booster 
cap cavity (Figure H-5) while the corresponding Mach number for the dual simultaneous firing 
case is very low (Mach number colors are blue in the booster cavity). The stagnation zone in the 
dual, simultaneous firing case will significantly reduce the convective heat transfer to the booster 
cap from the gas-phase despite the large values of temperatures that arise from compression of 
the gas (see Figure H-4(c) and (d)). Conversely, the convective heat transfer is enhanced for the 
single firing case due to the high temperature and Mach number values in the booster cavity 
(Figure H-3(c) and Figure H-5 (c)). The effect of the stagnation zone created by the dual, 
simultaneous firing on the particulate motion is described below. 

The concentration contours of the ZrC >2 particles are shown in Figure H-7 and Figure H-8. The 
complex shock structure that forms in the flame channels (see Figure H-l and Figure H-2) causes 
the particles to enter into the flame channel at an angle and slide along the inside wall. As 
discussed previously, the pressure in the booster cap region for the single firing case is allowed 
to relieve into the unused channel, whereas a stagnation zone is formed in the simultaneous, dual 
firing case eventually leading to backward flow. The effect of this difference on the particle- 
phase can be seen in the concentration contours at 150 ps in Figure H-7(d) and Figure H-8(d). 

For the single firing case, the particles are able to travel unimpeded into the booster cavity and a 
high amount of particles impact and collect on the booster cap. Conversely, the backward flow 
in the simultaneous, dual firing case impedes the forward progress of the particle phase and only 
trace amounts of particles are able to reach the booster cap potentially resulting in much reduced 
heat transfer and a less reliable burn through. Thus, the simulations presented here provide a 
credible mechanism to explain the failure mode in dual, simultaneous firings of the NSI 
initiators. 

H.3 Concluding Remarks 

The simulations performed on a 2-D planar configuration of the SS design indicate that the 
failure mode for simultaneous dual NSI firings may be caused by flow interactions between the 
flame channels. The shock waves from each initiator interact in the booster cap region resulting 
in a high pressure that prevents the gas and particulate velocity from rising in the booster cavity 
and in fact preventing the bulk of the particulate phase from impacting the booster cap. This 
reduces the heat transfer to the booster cap both because the particles do not reach the booster 
cap, and the heat transfer coefficient is reduced due to lack of convective effects. The current 
study will be continued to provide estimates of the heat transfer to the booster cap as well as the 
walls of the flow channels. In addition, 3-D simulations will be performed to further investigate 
the failure mode and ensure that the proposed theory suggested by the 2-D simulation continue to 
hold for the complete 3-D configuration. 
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(C) (d) 

Figure H-l. Pressure Distribution for Single Initiator Firing at 0, 50, 100, and 150 jus 
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(C) (d) 

Figure H-2. Pressure Distribution for Simultaneous Dual Firings at 0, 50, 100, and 150 jus 
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Figure H-3. Temperature Distribution for Single Initiator Firing at 0, 50, 100, and 150 ps 
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(a) (b) 



(c) (d) 

Figure H-4. Temperature Distribution for Simultaneous Dual Firings at 0, 50, 100, and 150 jus 
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(C) (d) 

Figure H-5. Mach Number Distribution for Single Firing at 0, 50, 100, and 150 jus 
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Figure H-6. Mach Number Distribution Simultaneous Dual Firing at 0, 50, 100, and 150 /.is 
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(C) (d) 

Figure H-7. Particulate Distribution for Single Firing at 0, 50, 100, and 150 jus 
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(a) (b) 


t = 0000150s 
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00 (d) 

Figure H-8. Particulate Distribution for Simultaneous Dual Firing at 0, 50, 100, and 150 /is 
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